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ABSTRACT 

In this paper we report the results of a 7 -ray study of IceCube’s extraterrestrial neu¬ 
trino candidates detected as track-like events. Using 70 months of Fermi -LAT obser¬ 
vations, a likelihood analysis of all 1 — 300 GeV photons within 5° of the track-like 
neutrino candidates’ origin was undertaken, to search for spatially coincident 7 -ray 
emission. One of IceCube’s HESE track events was found to be spatially coincident 
with a 7 -ray bright active galactic nucleus (AGN), PKS 0723-008. We find however, 
that the chance probability for Fermi -LAT detected AGN to be spatially coincident 
with a single HESE track-like event is high 37%). We therefore find no evidence 
of 7 -ray emission associated with the detection of IceCube’s HESE track-like neutrino 
candidates. Upper limits were calculated in the energy range of 1 — 300 GeV, assuming 
a point source origin for the neutrino events considered. The implications for the non¬ 
detection of 7 -ray emission from the source of the HESE track-like events are briefly 
discussed. 

The large time period analysed in our study did however, reveal two new 7 -ray 
point sources. With a flux of (1.28 =b 0.08) xlO -9 photons cm -2 s _1 , in the 1 — 300 
GeV energy range, and an associated TS value of 220.6, one of these new point sources 
is positionally coincident with the AGN PKS 1346-112. The other point source has 
a 1 — 300 GeV flux of (7.95 =b 1.23) xlO -10 photons cm -2 s -1 and an associated TS 
value of 92.4. This new point source is spatially coincident with the radio source NVSS 
J072534 + 021645 suggesting that it too is an AGN. 

Key words: radiation: non-thermal - galaxies: active - gamma rays: galaxies - 
neutrinos. 


1 INTRODUCTION 


In 2013 the IceCube Neutrino Observatory published obser¬ 
vational evidence for the existence of extraterrestrial high- 


energy neutrinos (IceCube Collabo ratio n (2013)). This evi¬ 
dence was obtained through a 4 HESE [^analysis of the first 
two years of data taken with the 79 and 86 string IceCube 
detector. The addition of a third year of data resulted in the 
evidence evolving to a 5 cr discovery of these extraterrestrial 


neutrinos (Aarsten et al. (2014a)). IceCube’s 3-year HESE 
analysis uncovered a total of 37 extraterrestrial neutrino 
candidates with deposited energy in the range of 30 TeV to 2 
PeV. The null hypothesis that all of these 37 events are from 
an atmospheric background was rejected at the 5.7cr level of 
confidence. Employing a series of likelihood-based tests, the 
IceCube collaboration found no significant evidence against 


* E-mail: anthony.brown@durham.ac.uk 

1 IceCube’s HESE analysis is its High Energy Starting Event 
analysis, aimed at selecting high energy events that have an in¬ 
teraction vertex within the IceCube detector. 


an isotropic distribution for the HESE neutrino candidates. 
With the discovery of extraterrestrial neutrinos, IceCube has 
open a new window through which to view the high-energy 
Universe. Due to their importance, these 37 candidate neu¬ 
trino events have been interpreted within a range of galactic 


and extra-galactic scenarios (e.g. 

. Razzaque 

(2013), Padovani 

V Resconi 

(2014), Dermer et al. 

|2014), Kohta, Inoue V Der- 

mer (2014 

), Krauss et al. (2014 

D- 


Neutrino emission from high-energy astrophysical 
sources capable of accelerating protons to relativistic ener¬ 
gies has long been predicted. The neutrinos arise from the 
decay of charged mesons created in the interactions of the 
relativistic protons with either ambient gas (p-p) or ambient 
radiation (p- 7 ). If the 20 TeV — 2 PeV extraterrestrial neu¬ 
trinos detected by IceCube were indeed created this way, 
then their existence implies the existence of astrophysical 
objects capable of accelerating protons to energies of at least 
0(10) to 0(100) PeV. The decay of neutral mesons created 
in the same interactions of the relativistic protons produce 
7 -rays. While the p-p interactions produce a smooth power 
law distribution for the resultant neutrinos and 7 -ray fluxes, 
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the p -7 interactions usually produce a peaked distribution 
for the resultant neutrinos and 7 -ray photons, with the peak 
and spread of these distributions strongly dependent upon 
the nature of the ‘target’ photon population (Miicke et al. 
(2000a,b)). Furthermore, the ratio of the typical neutrino 
and 7 -ray energy, E 1 /E„ , also depends upon the character¬ 
istics of the target photons, with a ratio of ~ 1.2 expected 
for a synchrotron radiation target and ~ 1.1 for a thermal 
radiation target distribution. Regardless of whether the in¬ 
teraction is synchrotron p -7 or thermal p- 7 , the ~ 10 — 100 s 
TeV 7 -rays associated with the ~ 30 — 250 TeV neutrino 
candidates considered in this study are believed to quickly 
cascade down to E 1 < 100 GeV energies due to the photon 
opacity of the emission region ( |Mannhe im (1993)). Assum¬ 
ing that the opacity of the neutrino/ 7 -ray emission region 
is large enough to facilitate the transformation of the initial 
TeV photons down to GeV energies while not being so large 
as to absorb the 7 -ray flux entirely, a search for the neutrino 
sources can be made with the Fermi- LAT detector. 

The LAT detector onboard the Fermi satellite, de¬ 


scribed in detail by jAtwood et al. (2009), is a pair-conversion 


telescope, sensitive to a photon energy range from below 20 
MeV to above 300 GeV. With a large field of view, ~ 2.4 
sr, improved angular resolution, 68 % containment angle of 
~ 0.8° at 1 GeV, and large effective area, ~ 8000 cm 2 
on axis for 10 GeV photons, Fermi- LAT provides an or¬ 
der of magnitude improvement in performance compared to 
its EGRET predecessor. Since 2008 August 4, the vast ma¬ 
jority of data taken by Fermi has been performed in all- 
sky-survey mode, whereby the Fermi- LAT detector points 
away from the Earth and rocks north and south of its orbital 
plane, on subsequent orbits. This rocking motion, coupled 
with Ferraz-LAT’s large effective area, allows Fermi to scan 
the entire 7 -ray sky every two orbits, or approximately every 
three hours. This observational characteristic affords us, for 
the first time, continuous monitoring of the high-energy 7 - 
ray sky and thus has important implications for the study of 
transient events. In addition to this, coupling Fermi- LAT’s 
continual scanning of the sky with a long mission lifetime 
allows us to construct a deep exposure of the 7 -ray sky. 

IceCube classifies events by the pattern of the 
Cherenkov light seen in the detector array; track events 
are those where the presence of a high-energy muon has 
produced a distinctive visible track of light while shower 
events are those where the Cherenkov light, from the parti¬ 
cle shower initiated when the neutrino interacts, is observed 
only in aggregate. The angular resolution for events con¬ 
taining visible muon tracks is much better 1 °) than for 
those that do not 15° (Aarsten et al. ( |2014a| )). For equal 
neutrino fluxes of all flavours, events with tracks make up 


only 20% of interactions (Aarsten et al. (2014a)). Eight of 
the 37 HESE neutrino candidate events were classified as 
track events, 28 as shower events and one event was formed 
by two coincident cosmic-ray muon events. Although there 
are fewer track events, the angular resolution of the track 
events is comparable with that of the Fermi- LAT detec¬ 
tor for E 1 > 1 GeV photons which facilitates associating 
the track-like events of IceCubes HESE analysis with known 
Fermi- LAT detected F# > 1 GeV 7 -ray sources. 

An important product of Ferraz-LAT’s continual scan¬ 
ning of the sky is the 7 -ray catalogues that it periodically 
produces such as the one and two year point source cat¬ 


alogues (1FGL (Abdo et al. (2010a)) and 2FGL (Nolan 
|et al.| p(H2ll)), or the second Active Galaxies catalogues 
(2LAC ( jAckermann et ah ( 2011| )). The Ferraz-LAT col¬ 
laboration recently released their first high-energy catalog 
( 1 FHL), which lists all F# >10 GeV sources from the first 


three years of Fermi- LAT observations (Ackermann et al. 


( 2011 }) 


Padovani V Resconi (2014) looked for astrophysical 


counterparts to the E u > 60 TeV HESE neutrino candi¬ 
dates in three 7 -ray catalogues including the 1FHL cata¬ 
logue. Padovani V Resconi placed their highest priority on 
the TeVCat catalogue which reaches the TeV regime. Un¬ 
surprisingly, given that the neutrino sample of their study 
was dominated by shower events (83%) with an average me¬ 
dian angular error of ~ 11.4°, they found a large amount 
of spatial correlation between the catalogues and their neu¬ 
trino candidates. To therefore determine the plausibility of 
the association between the TeV emitting coincident objects, 
they had identified, and the neutrino candidates, Padovani 
V Resconi compared the 7 -ray spectral energy density of the 
object with the neutrino flux per event at the specific en¬ 
ergy of the relevant neutrino and assessed whether a simple 
extrapolation plausibly connected the 7 -rays with the neu¬ 
trino. Around 10 % of the TeVCat counterparts passed this 
energy diagnostic. For the counterparts which were listed 
only in the 1FHL catalogue, Padovani V Resconi decided 
there was too big a gap between the source’s highest de¬ 
tected 7 -ray energy and the neutrino energy to make a sen¬ 
sible link. Padovani V Resconi found no spatial correlation 
between the 1FHL catalogue and the median angular error 
radius of the three track-like events in their neutrino sample, 
although they did note that the 2 FGL source PKS 0723-008 
is correlated with neutrino candidate #5 and that the 1FHL 
source 4C+41.11 and the 2 FGL source PKS 0823-321 were 
just outside the median angular errors of candidates #13 
and #3 respectively. 

In what could be considered the inverse cross¬ 
correlation, the IceCube collaboration have recently per¬ 
formed a likelihood search for cumulative neutrino emission 


from a large sample of of 2 LAC sources (Glusenkamp et al. 


(2015)). The neutrino sample used was a muon-track sam¬ 
ple from the years 2009-2011 containing over 300,000 events 
with a median angular resolution of 1° at 1 TeV and 0.5° 
at 100 TeV. In the southern hemisphere the sample is dom¬ 
inated by muons from air showers and in the upgoing re¬ 
gion the majority of events come from atmospheric muon 
neutrinos. This was accounted for in the likelihood analysis 
through the background probability distribution function. 
The study found no significant correlation and concluded 
that 2 LAC blazars are not responsible for the majority of 
the diffuse astrophysical neutrinos reported in [Aarsten et al.| 
(2014a). 

Thus, there is currently, no clear association between 
the IceCube extraterrestrial neutrino candidates and 7 -ray 
sources and as Padovani V Resconi remarked, this strongly 
motivates more TeV photon observations. If it is the case, 
as mentioned above, that the ~ 10 — 100s TeV 7 -rays 
quickly cascade down to F# < 100 GeV energies then 
the 7 -ray counterparts of the neutrino sources may be re¬ 
vealed through a close analysis of Ferraz-LAT observations. 
Here, rather than cross-correlating against precompiled cat¬ 
alogues, we use the deep exposure of the first 70 months 
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of Fermi -LAT operation to search for GeV counterparts to 
the track events of IceCube’s 3-year HESE analysis. Such 
an analysis, while sensitive to catalogued sources, also al¬ 
lows us to search for faint 7 -ray sources not present in the 
existing Fermi catalogues and, in the absence of evidence 
for a source, determine upper limits for the associated 7 - 
ray flux. Also the 70 month Fermi -LAT data set contains 
the 3 year time period covered by the IceCube HESE search 
which allows our analysis to be sensitive to the possibility 
that the neutrinos were produced during the flare of a nor¬ 
mally faint source. All neutrino candidates with a track-like 
topology are investigated, with the exception of event #28, 
which displays evidence of a cosmic ray induced air shower in 
the IceTop detector accompanying the neutrino candidate, 
thus suggesting it is part of the atmospheric background. 
Each of the remaining neutrino candidates are considered 
independently. In Section 2, we describe the Ferrm-LAT ob¬ 
servations and data analysis routines employed during our 
study. Section 3 touches on the GeV properties of the new 
7 -ray point sources discovered during our study, with Sec¬ 
tion 4 briefly discussing possible interpretations of the non¬ 
detection of GeV emission spatially co-incident with the ori¬ 
gin of the neutrino candidates considered. 


2 FERMI-LAT OBSERVATIONS AND DATA 
ANALYSIS 

The data used in this study comprise all Fermi -LAT event 
and spacecraft data taken during the first 70 months of 
Fermi -LAT operation, from 2008 August 4 to 2014 June 
4, which equates to a Mission Elapsed Time (MET) inter- 
val of 239557417 to 423589417. All SOURCE 7 -ray event^ 
in the 1 < E 1 < 300 GeV energy range, within a 5° radius 
of interest (Rol) centered on the right ascension and decli¬ 
nation of the origin of HESE track-like events were consid¬ 
ered In accordance with the PASS7_REP criteria, a zenith 
cut of 100 ° was applied to the data to remove any cos¬ 
mic ray induced 7 -rays from the Earth’s atmosphere. The 
good time intervals were generated by applying a filter ex¬ 
pression of “(data_qual==1) &;& (lat_config==1) &&; 
ABS(rock_ANGLe)< 52” to the dats[^] A summary of the 
analysis criteria is given in Table ^ 

Throughout this analysis, version v9r33p0 of the 
Fermi Science Tools was used in conjunction with the 
p7rep_SOURCE_v15 instrument response functions (IRFs). 
During the likelihood analysis, a model file consisting of both 
point and diffuse sources of 7 -rays was employed. For the dif¬ 
fuse emission, the most recent Galactic, gll_iem_v05_revl.fit, 
and extragalactic, iso_source_v05.txt, diffuse models were 
used. In addition to these diffuse models, all point sources 
within 10 ° of the candidate neutrino position (aj 2 ooo, 


2 SOURCE events have an event class of 2 in the pass 7_REP data, 


al. 

(2012a) 

3 r 

l’he rign 


events were taken from Table 1 of|Aarsten et al.|( 2014a). 

4 See http: //f ermi .gsf c .nasa.gov7ssc7d.ata/analysis/ 

documentation/Cicerone/Cicerone_Data_Exploration/ 
Data_preparation.html for details on LAT data selection. 


Table 1. Summary of the criteria utilised in this analysis. 


Science Tools version 
IRF 

Event class 
Photon Energies 
Rol 

Zenith angle cut 
Rocking angle cut 
LAT config/Data quality 
Galactic diffuse model 
Isotropic diffuse model 


v9r33p0 

p7rep_source_v15 
source, Reprocessed Pass 7 
1 - 300 GeV 
5° 

< 100 ° 

< 52° 

gll_iem_v05_revl .fit 
iso_source_v05.txt 


$J 2 ooo) were considered in the likelihood model. The po¬ 
sition and spectral shape for these point sources were taken 
from the 2FGL. The normalisations and spectral indices of 
all point sources in the model file were left free to vary, with 
the exception of those within the 5° to 10° annulus from 
the candidate neutrino event (<aj 2 ooo, < 5 j 2 ooo) position which 
were fixed to their 2FGL catalogue values. The normalisa¬ 
tion factor of the extragalactic diffuse emission was left free 
to vary, while the Galactic diffuse template was multiplied 
by a power law in energy, the normalisation of which was 
also left free to vary. The model map immediately allows us 
to determine if there are any catalogued sources spatially 
coincident with the neutrino candidates and, as noted by 


Padovani V Resconi (2014), there is one such coincidence 


with the 2FGL AGN located within median angular error of 
neutrino candidate #5. The significance of this spatial cor¬ 
relation, along with the observed 7 -ray properties over the 
70-months, are investigated in detail in §2.1. 

An initial binned likelihood analysis was conducted to 
confirm the accuracy of the ‘diffuse + 2FGL point source’ 
model. The best-fit model from this initial likelihood fit was 
used, in conjunction with the Fermi tool GTMODEL, to con¬ 
struct a ‘model map’ of a 7°x7° region centered on the 
J 2000 1 ^J 2 ooo) of each candidate neutrino track event. This 
model map was then compared with the 70-month sky map 
of the 1 — 300 GeV 7 -rays to create a residuals count majJj 
Any positive excess of events in the residuals map indicates 
the presence of a new source of 7 -ray^] If new sources are 
present in the data, and not properly accounted for within 
the model file utilised during the likelihood analysis, they 
can artihcally increase the significance of the 7 -ray flux from 
the neutrino candidates. The 7°x7° sky map, model map 
and residuals maps for the 7 separate HESE track events 
can be seen in Figures 1 to 7. All maps have been smoothed 
with a 1° Gaussian. The colour scales for the sky and model 
maps are in units of 7 -ray counts, while the residuals maps 
are in units of percentage. 

For the majority of the candidate neutrino events, the 
residuals maps do not reveal any obvious excess in 7 -ray 
emission coincident with the neutrino’s point of origin. In¬ 
deed, the residuals maps for events #3, # 8 , #13, #18 and 


5 The residuals count map was constructed via (sky map—model 
map)/model map (e.g. see " 


Brown et al. 


(2014)) 


6 The use of a residuals map" to look for new 7 -ray sources is 
independent and complementary to the alternative approach of 
using the GTTSMAP Fermi tool (see eg. |Brown &; Adams ( 2 012| ), 


Macias et al. 


( 2012 ) & 


Brown 


(2014)). 
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188 1.1 2.2 3.3 4.3 5.4 6.4 7.5 8.5 9.6 1 .1 2 2.9 3.8 4.7 5.6 6.5 7.4 8.3 9.2 1 95 -0.65 -0.34 -0.034 0.27 0.58 0.88 1.2 1.5 1.8 2 

Figure 1 . Event 3: Left: 7°x7° sky map of the 1 — 300 GeV flux for the entire 5.8 year data set. Middle: 7°x7° 
model map produced with the Fermi tool GTMODEL, utilising the best-fit model from the 5.8 year binned likelihood 
analysis. Right: 7° x7° residuals map of the 1 — 300 GeV events. The residuals map is produced by (sky map—model 
map)/model map). All maps have been smoothed with a 1° Gaussian. The cross at the centre of each map indicates 
the (ckJ 20005 <5j200o) of the neutrino candidate, while the circle present in each map indicates the median error 
radius for the neutrino candidate as shown in Table 2. The colour scales for the sky and model maps are in units of 
7 -ray counts, while the residuals maps are in units of percentage.There is no evidence for new 7 -ray sources beyond 
the diffuse + 2FGL point sources model used. 



3 1.9 3.8 5.7 7.6 9.5 11 13 15 17 1 95 3.2 5.4 7.6 9.8 12 14 16 19 21 2 1 -0.52 -0.035 0.45 0.93 1.4 1.9 2.4 2.9 3.4 3 


Figure 2. 7°x7° sky, model and residuals maps for Event 5. There is an excess in the residuals map consistent 
with a point source at (<aj 2 oocb ^J 200 o) — (111.419, 2.221). 



3 1.1 2.1 3.2 4.2 5.3 6.4 7.4 8.5 9.5 1 78 1.8 2.7 3.7 4.7 5.7 6.7 7.7 8.6 9.6 1 1 -0.54 -0.078 0.39 0.85 1.3 1.8 2.2 2.7 3.2 3 


Figure 3. 7°x7° sky, model and residuals maps for Event 8 . There is no evidence for new 7 -ray sources beyond 
the diffuse + 2FGL point sources model used. 
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61 4.4 8.2 12 16 20 23 27 31 35 3 .6 7.1 12 16 21 25 30 34 39 44 4 86 -0.43 -0.0083 0.42 0.84 1 


Figure 4. 7°x7° sky, model and residuals maps for Event 13. There is no evidence for new 7 -ray sources beyond 
the diffuse + 2FGL point sources model used. 



3 2.9 5.9 8.8 12 15 18 21 24 26 2 53 3.7 6.8 9.9 13 16 19 22 26 29 3 1 -0.49 0.029 0.55 1.1 1.6 2.1 2.6 3.1 3.6 4 


Figure 5. 7°x7° sky, model and residuals maps for Event 18. There is no evidence for new 7 -ray sources beyond 
the diffuse + 2FGL point sources model used. 


#37, show no evidence of additional point sources beyond 
the 2FGL point sources described in the model file. However, 
the residuals maps for events #5 and #23, seen in Figures 2 
and 6 respectively, show a clear point source excess. These 
excesses indicate the presence of a new 7 -ray point source, 
which is not in the 2 FGL catalogue, in the field of view. 
These new 7 -ray sources will be discussed in detail in §3. 

To search for faint 7 -ray sources associated with the 
candidate neutrino’s origin, that may not be apparent in the 
residual sky map, a further BINNED analysis was conducted, 
with an additional point source included in the model file at 
the position of the candidate neutrino event. The spectral 
form of these additional point sources was described by a 
power law spectrum, dN/dE = A x(E/E a )~ r , with both A 
and T left free to vary. The (<aj 2 ooo, $J 2 ooo) for this ‘neu¬ 
trino’ point source were fixed to the values published by 
the IceCube Collaboration (IceCube Collaboration (2013) 
& Aarsten et al. (2014a)). 

All ‘neutrino’ point sources added to the model file for 
this additional 70-month BINNED analysis had a test statis- 
ticQof 0. As such, no evidence for 1 — 300 GeV 7 -ray emission 


7 The test statistic, TS, is defined as twice the difference between 
the log-likelihood of two different models, TS = 2 [logL — logLo], 
where L and Lq are defined as the likelihoods when the source is 
included or not respectively ( [Mattox et a l. (1996|). 


spatially coincident with the (<aj 2 ooo, < 5 j 2 ooo) of the HESE 
track events was found. Upper limits, at the 95% confidence 
level, were therefore calculated for the 7 -ray emission, as¬ 
suming a point source for the candidate neutrino events. 
These upper limits are summarised in Table 2. 


2.1 Icecube neutrino event #5 


In addition to a new 7 -ray source within the field of view 
studied, the HESE event #5 has a known 7 -ray bright AGN 
within the la uncertainty of the event’s position of origin. 
Located (1.04±0.04)° from event #5 is PKS 0723-008. PKS 
0723-008 is a BL Lac object, present in both the 1 year 
and 2 year Ferrm-LAT AGN catalogues (1LAC and 2LAC; 
Abdo et al. (2010b), Ackermann et al. (2011)). In our 70- 
month binned analysis, PKS 0723-008 is detected with a 
TS of 176.8, a spectral index T m 2.07 =b 0.12 and a flux of 
(9.35 ± 1.20) x 10 - 10 ph cm -2 s _1 . Both PKS 0723-008’s 70- 
month flux and spectral index are consistent with the values 
published in the 2 LAC catalogue. 

To quantify the significance of the spatial coincidence 
of 2LAC AGN with a single track-like neutrino event, the 
positions of all 2 LAC sources were scrambled in aj 2 ooo and 
$J 2 ooo, and the new positions were cross-correlated with the 
<aj 2 ooo and (U 2000 of all neutrino candidates with E v ^ 60 
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Figure 6. 7°x7° sky, model and residuals maps for Event 23. There is a excess in the residuals map consistent 
with a point source at (aj 2 oocb ^J 200 o) — (207.409,-11.553). 



3 1.7 3.5 5.2 7 8.7 10 12 14 16 1 53 2.5 4.5 6.4 8.4 10 12 14 16 18 2 1 -0.38 0.25 0.88 1.5 2.1 2.8 3.4 4 4.6 5 


Figure 7. 7°x7° sky, model and residuals maps for Event 37. There is no evidence for new 7 -ray sources beyond 
the diffuse + 2FGL point sources model used. 


Te\|^] considered in this study. Repeated 10,000 times, this 
approach revealed that the chance probability of having a 
2LAC AGN positionally coincident with one of the E u ^ 60 
TeV track-like neutrino candidates to be ~ 37%. This large 
probability suggests that the spatial coincidence between 
event #5 and PKS 0723-008 is a chance event. 

Nonetheless, to investigate the temporal characteris¬ 
tics of PKS 0723-008’s 1 — 300 GeV 7 -ray flux around the 
time the #5 neutrino candidate was detected, the data were 
binned into 50 day periods, with an UNBINNED GTLIKE rou¬ 
tine applied to each bin separately. Only temporal bins 
with test statistic TS > 10 were considered. The resultant 
lightcurve is shown in Figure [ 8 ] While the lightcurve does re¬ 
veal periods of enhanced 7 -ray flux from PKS 0723-008 when 
compared to the 70-month average, there is no obvious flare 
event associated with the detection of the candidate neu¬ 
trino event. Indeed, during the 50-day period in which the 
neutrino candidate was detected by IceCube, no significant 
7 -ray flux was observed. 

Phenomenological studies by Fermi have found the ma¬ 
jority of TeV bright AGN to have T < 2 (Abdo et al. (2009)). 
Since the 70-month averaged spectrum of PKS 0723-008 is 


8 As shown in §4, the background neutrino events considered in 
this study are expected to concentrate at lower energies, E u < 60 
TeV. 


marginally softer than this critical value, we searched for 
E 1 > 100 GeV photons spatially coincident with PKS 0723- 
008. A closer inspection of the individual photon events re¬ 
veals the presence of a 115 GeV ultraclean photon within 
0°.l of PKS 0723-008. Utilising the ‘diffuse + 2FGL + new 
source’ model, with all model parameters frozen to their 70 
month best fit parameters, the Fermi tool GTSRCPROB calcu¬ 
lated the probability that the VHE photon originated from 
PKS 0723-008 to be 0.97806. Assuming Gaussian errors, this 
equates to < 2 . 5 cr significance of the 115 GeV photon orig¬ 
inating from PKS 0723-008. There is therefore no evidence 
of PKS 0723-008 being a source of VHE photons. 


3 NEW 7 - RAY POINT SOURCES 

The 70-month long exposure analysed in this study revealed 
two new 7 -ray point sources in the field of view of IceCube’s 
HESE events #5 and #23. To investigate the global 7 -ray 
properties of these new point sources, an additional BINNED 
analysis was performed for both of these events, with the 
model file modified to include an additional point source, 
described by a power law and spatially coincident with the 
excesses in residuals map plots of Figures 2 & 6 . 

For the point source in the field of view of the #5 event, 
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Table 2. Summary of upper limits for the 1 — 300 GeV 7 -ray flux spatially coincident with the (<aj 2000 , ^J 200 o) of 
the HESE track events. 


IceCube event 

Deposited Energy 
(TeV) 

“J2000 

(deg) 

1/2000 

(deg) 

Angular error 
(deg) 

Northern/Southern 

Hemisphere 

7 -ray flux upper limit 
(ph cm -2 s -1 ) 

#3 

78.7lg°f 8 

127.9 

-31.2 

7 1.4 

South 

1.89 X 10 “ 10 

#5 

71 4+ 9 - 0 
‘ 1 - 9.0 

110.6 

-0.4 

7 1.2 

South 

9.09 x 10“ n 

#8 

32.6R°; 8 

182.4 

- 21.2 

^ 1.3 

South 

8.64 x 10“ n 

#13 

253+2® 

67.9 

40.3 

7 1.2 

North 

1.99 x 10 “ 10 

#18 

3i.5±^:i 

345.6 

-24.8 

< 1.3 

South 

1.05 x 10 “ 10 

#23 

82.2+g'® 

208.7 

-13.2 

< 1.9 

South 

1.32 x 10 “ 10 

#37 

30.81!;! 

167.3 

20.7 

7 1.2 

North 

5.18 x 10“ n 



Figure 8. The lightcurve of the 1 — 300 GeV flux from PKS 0723- 
008 binned in 50 day periods. The horizontal dashed line shows 
the flux from the 70-month binned analysis, while the vertical 
solid line shows the time the #5 neutrino candidate was detected. 
While there are periods of enhanced 7 -ray flux from PKS 0723- 
008 when compared to the 70-month average, there is no obvious 
flare event associated with the detection of the neutrino event. 


the resultant best fit was found to be: 


dN 

dE 


= (1.01 ±0.24) x 10 -14 


E 

6938 MeV 


-2.441±0.189 


photons cm 2 s 1 MeV 1 ( 1 ) 


which equates to an integrated flux of 
T± 7 >i GeV = (7.95 ± 1.23) x 10 -10 photons cm“ 2 s _1 (2) 


taking only statistical errors into account. The test statistic 
value of the best-fit power law was TS = 92.4, equating to 
a significance of ~ 9.6cr. 

For the point source in the field of view of the #23 
event, the resultant best fit was found to be: 


dN 

dE 


(1.761 ±0.01) x 10“ n ( 


E 


\ —2.752 i 0.003 


473 MeV 


photons cm _ 2 s _ 1 MeV _1 (3) 

which equates to an integrated flux of 

Fe i >1 GeV = (1.28 ±0.08) x 10 9 photons cm 2 s 1 (4) 

again taking only statistical errors into account. The TS 
value of the best-fit power law was TS — 220.6, equating to 
a significance of ~ 14.9cr. 

The association of these new 7 -ray sources with known 
astrophysical sources requires a more precise localisation 


of the 7 -ray emission. This was done with the Fermi tool 
GTFINDSRC, in conjunction with an unbinned exposure map. 
Using the best-fit model from the 70-month analysis, the 
GTFINDSRC derived (<aj 2 ooo, ± 72000 ), along with the la error 
radius, for the two new point sources, are given in Table [3] 
The new 7 -ray point source within event #23 ’s field 
of view is spatially coincident with the flat spectrum ra¬ 
dio quasar (FSRQ) PKS 1346-11^] For the new 7 -ray point 
source within event #5’s field of view, while there is no coun¬ 
terpart within GTFlNDSRC’s the la error radius, the radio 
source NVSS J072534 ± 021645 does lie within the 2a error 
radius of the new source. This spatial coincidence suggests 
that this new 7 -ray point source is a blazar. 


4 DISCUSSION 

Of all extraterrestrial neutrino candidates considered in this 
study, only event #5 has a 7 -ray point source within the 
median angular error of the track reconstruction, with the 
BL Lac object PKS 0723-008 located ~ 1° from #5’s origin. 
However, Monte Carlo simulations revealed the probability 
that a 2LAC source is positionally coincident with event #5 
to be ~ 37%, suggesting that this positional coincidence is 
merely a chance event. 

In addition to the high probability that a 2LAC source 
is positionally coincident with a single track-like neutrino 
event, there are no distinguishing characteristics in the 
1 — 300 GeV 7 -ray properties, when compared to the other 
BL Lac objects in the 2LAC catalogue, to suggest PKS 
0723-008 is a source of neutrinos. For example, while the 
50-day binned lightcurve shown in Figure 8 , reveals peri¬ 
ods of increased 7 -ray flux (2 — 3) times greater than the 
70-month average, these flare periods do not coincide with 
the detection of neutrino event #5. Indeed, during the 50 
day period in which the neutrino candidate was observed, 
no significant 1 — 300 GeV 7 -ray emission was observed from 
PKS 0723-008. Likewise, while PKS 0723-008 does appear to 
have a fairly hard spectrum, there is no evidence that PKS 
0723-008 is a source of VHE photons. So in summary the 
combination of the insignificant spatial coincidence and the 


9 We note that the recently published third Fermi source cata¬ 
logue has found PKS 1346-112 to be 7 -ray bright from 4 years of 
data (Acero et al. (2015)). 
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Table 3. Summary of (<aj 2 oocb ^J 200 o ) 5 flux and spectral indices for the new 7 -ray point sources detected in the 
field of view for event #5 and #23. The new point source in the field of view of #5 is positionally coincident with 
the radio source NVSS J072534±021645 while the point source in the field of view of #23 is positionally coincident 
with the AGN PKS 1346-112. 


Event # 

1 - 300 GeV Flux 
(ph cm -2 s * 1 ) 

Spectral Index 

TS 

«J2000 

(deg) 

$j 2000 
(deg) 

lcr error radius 
(deg) 

Associated source 

5 

(7.95 ± 1.23) xlCT 10 * * * * * 

2.441 ±0.189 

92.4 

111.436 

2.23478 

0.039 

NVSS J072534 ± 021645 

23 

(1.28 ±0.08) xlO -9 

2.752 ± 0.003 

220.6 

207.411 

-11.5502 

0.031 

PKS 1346-112 


the observed 7 -ray properties suggest that PKS 0723-008 is 
not the source of the #5 neutrino candidate. Therefore, this 
study has found the origin of the HESE track events to be 
7 -ray dark in the 1 — 300 GeV energy range. 

A possible explanation for the non-detection of 1 — 300 
GeV 7 -ray emission spatially coincident with the observed 
neutrino events is that the diffuse neutrino flux observed by 
IceCube is the aggregate of a very large number of faint neu¬ 
trino sources which, correspondingly, are faint 7 -ray sources. 
IceCubes own point source limits (Aarsten et al. (2014b)) 
mean that a diffuse astrophysical flux at the leve l of 10~ 8 
GeV cm -2 s -1 sr -1 , as reported in (Aarsten et al. (2014a)), 
must be distributed over at least 0(50) sources. In order to 
be able to relate the neutrino flux from an individual source 
to a corresponding photon flux at GeV energies requires de¬ 
tailed modelling which is beyond the scope of this paper. 


quickly be processed to energies well below the photon en¬ 
ergy range considered in this studw Observations with the 
future Cherenkov Telescope Arra)fjor HAWC (Acharya et 


However we note that typical models (Saba et al. (2013), 
Doert et al. (2012), Halzen V Zas (1997)) relating a source’s 

photon and neutrino bolometric luminosities, would suggest 
that a photon flux, at the level of the upper limits presented 
in Table 2, is associated with an E ~ 2 neutrino flux at least an 
order of magnitude below the current IceCube point source 
limits. Thus, if the TeV photons associated with the neutral 
mesons are able to quickly cascade down to the GeV energies 
considered in this dedicated analysis of Fermi -LAT obser¬ 
vations and if the typical relationships between the neutrino 
and photon bolometric luminosities hold, then our photon 
limits suggest that the diffuse neutrino flux observed by Ice- 
Cube is drawn from the aggregate of at least 0(500) sources. 

Alternative explanations could be that the assumption 
of efficient cascade of TeV photons to GeV energies is not 
valid, or that the GeV emission is absorbed after leaving 
the neutrino source or that the neutrino candidates we have 
studied are in fact residual atmospheric background events. 
We briefly touch on each of these possibilities below. 

(i) The ratio of the neutrino and 7 -ray energies pro¬ 

duced through the decay of 7r°, 7r + and mesons is close 
to unity, such that the 30 — 250 TeV neutrinos considered 

in this study should be accompanied by 7 -rays of a similiar 

energy. These TeV photons are believed to quickly cascade 

down to <100 GeV energies due to the photon opacity of 
the emission region (Mannheim (1993)). In the context of 
the emission region’s 7 -ray opacity, the 7 -ray non-detection 

of the neutrino sources of our study can be interpreted as 

being due to the opacity being either too small or too large. 

If the opacity is too low, the TeV 7 -rays are able to freely 
escape and thus the 7 -ray flux is not able to cascade down 
to the photon energy range considered in this study. Like¬ 
wise, if the opacity is too high, the TeV 7 -ray flux would 


al. (2013); Abeysekara et al. (2013)) would shed light onto 


the former possibility, while detailed multi-wavelength stud¬ 
ies would be needed to investigate the latter possibility. 

(ii) In parallel to the effects of internal absorption, 
the lack of 7 -ray emission associated with the HESE track 
events could be the result of photon absorption external 
to the neutrino source. This attentuation primarily oc¬ 
curs through the process of photon-photon pair production 
(77 -7 e + e _ ). For 7 -rays, the photon field primarily respon¬ 
sible for this photon-photon attentuation is the Extragalac- 
tic Background Light (EBL; e.g |Fazio V Stecker| (1970|). A 
study undertaken by the HESS collaboration found that the 
intensity of the EBL, at near-infrared/optical wavelengths, 
is close to the lower limit given by the integrated light of re¬ 
solved galaxies ( Aharonian et al.| ( 2006| ). More recently, the 
Ferrm-LAT collaboration extended this study using 1 — 500 
GeV data from the first 46 months of Ferrm-LAT operation, 


and drew a similiar conclusion (Ackermann et al. ( 2012 c)). 
Importantly, the Ferrm-LAT collaboration also defined a 
critical photon energy, E cr i t , below which ^ 5% of the 7 -ray 
photons are absorbed by the EBL. The lowest E cr n they de¬ 
rived was ~ 30 GeV for 7 -ray sources in the redshift range 
0.5 < z < 1 . 6 . As such, while invoking 77 absorption by the 
EBL as the primary reason for the 7 -ray non-detection can 
account for the lack of E 1 > 30 GeV flux, it cannot solely 
explain the lack of 1 — 30 GeV flux. 

(iii) Finally, we note that IceCube’s HESE analysis 
claimed the discovery of extraterrestrial neutrinos at a ~ 
5.6cr confidence level, with 37 candidate neutrino events de¬ 
tected relative to an expected background of 8.4±4.2 cosmic 
ray muon events and 6.61^’6 atmospheric neutrinos. As such, 
anywhere between a quarter and two-thirds of the neutrino 
candidates detected are expected to be background events 
from cosmic rays interacting with the Earth’s atmosphere. 
These background events are expected to concentrate at 
lower energies (see Figure 2 of Aarsten et al. (2014a)). In 
Padovani V Resconi (2014), the residual atmospheric back¬ 


ground contamination was reduced by only considering the 
E u ^ 60 TeV neutrino candidates in their correlation stud¬ 
ies. Applying an E u ^ 60 TeV cut to the track events of 
our study removes 3 neutrino candidates, events # 8 , #18 


10 It is worth highlighting that the upper limits of Table 2 are 
comparable to the 50 hour (20 GeV - 70 TeV) integrated sensi¬ 


tivity of CTA (Bernlohr et al. (2013)). 
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and #37. A closer inspection of the background atmospheric 
rate^j expected in the 30 — 40 TeV energy range and dec¬ 
lination bands of these removed events, as shown by sup¬ 
plemental Figure 5 of Aarsten et al. (2014a), reveals an ex¬ 
pected background of 1-2 events in the 10 ° declination band 
around events #8 and #18 and around 1 event in the 10 ° 
declination band around event # 37, making it rather likely 
that these neutrino candidates have an atmospheric origin. 

However, inspecting the expected background rates for 
the E u ^ 60 TeV neutrino candidates in their respective 
energy and declination band reveals a smaller number of ex¬ 
pected background atmospheric events, ranging from < 0.1 
to 0.25. In particular for event #13, which has a deposited 
energy of 253^22 TeV and a northern hemisphere origin, 
~ 0.06 background events are expected in its declination 
band during the 988 day IceCube livetime of the 3 year 


HESE analysis (see supplemental Figure 5 of Aarsten et al. 


(2014a)). This means that, while a contribution from the 
residual atmospheric background contamination could ex¬ 
plain some of the non-detections of possible 7 -ray counter¬ 
parts to the HESE track events, particularly at lower ener¬ 
gies, it cannot explain all of them. 


to transient events, the survey nature of HAWC observations 
would be required. These observations are strongly encour¬ 
aged. 
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5 CONCLUSIONS 

The angular resolution of a E 1 « 1 GeV photon detected by 
the Fermi -LAT detector is comparable to the median angu¬ 
lar error of IceCube’s HESE neutrino candidates detected as 
track-type events. As such, associating the track-like events 
of IceCube’s HESE analysis with known Fermi -LAT de¬ 
tected U 7 > 1 GeV sources is simplified greatly when com¬ 
pared to trying to find counterparts to IceCube’s shower-like 
events. Motivated by this, we analysed 70 months of Fermi- 
LAT observations to search for 1 — 300 GeV emission spa¬ 
tially coincident with the origin of IceCube’s HESE track 
events. With the exception of event #5, no E 1 > 1 GeV 
emission spatially coincident with the origin of the neutrino 
events was found. For these events, 95% confidence level up¬ 
per limits were calculated. 

Event #5 is the only HESE track-like event that has 
a 7 -ray point source within the median error of its origin, 
with PKS 0723-008 being separated from the origin by ~ 1 °. 
A Monte Carlo study found that the probability of such 
a spatial coincidence occuring by chance is relatively high, 
with it occuring in 37% of the cases we simulated. As such, 
an insignificant spatial coincidence, combined with the ob¬ 
served 7 -ray properties, suggests that PKS 0723-008 is not 
the source of the #5 neutrino candidate. 

The study undertaken here has been conducted on the 
assumption that the intrinsic TeV photons associated with 
the neutral mesons are able to quickly cascade down to GeV 
energies without being absorbed entirely. We note that, if 
this assumption is incorrect, and the TeV photons are able 
to freely escape from the source, observations with ground- 
based arrays such as HESS, VERITAS and CTA are required 
to search for the 7 -ray counterparts to IceCube’s neutrino 
events. However, if the production of these neutrinos is due 

11 Here we refer to the atmospheric background as being atmo¬ 
spheric muons and atmospheric neutrinos. 
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